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1 Summary 

The objective of this work has been to produce methodologies for high speed jet noise 
reduction based on natural mechanisms and enhanced feedback control to affect fre- 
quencies and structures in a prescribed manner. In this effort the two-point hot wire 
measurements obtained in the Langley jet facility by Ukeiley were used in conjuction 
with linear stochastic estimation (LSE) [1] to implement the LSE component of the 
complementary technique (v. Bonnet, Cole, Delville, Glauser and Ukeiley 1994, Ex- 
periments in Fluids). This method combines the Proper Orthogonal Decomposition 
(POD) [4] and LSE to provide an experimental low dimensional time depen- 
dent description of the flow field. From such a description it should be possible to 
identify short time high strain rate events in the jet which contribute to the noise (v. 
J.M. Seiner 1997, AIAA Fluid Dynamics Conference in Snowmass Village, CO). The 
main task completed for this effort is summarized below’: 

LSE experiments were performed at the downstream locations where the two point 
hot wire measurements have been obtained by Ukeiley. These experiments involved 
sampling simultaneously hot wire signals from a relatively course spatial grid in r 
and 6. From this simultaneous data, coupled with the two-point measurements of 
Ukeiley via the LSE components of the complementary technique, an experimental 
low dimensional description of the jet at 4, 5, 6, 7 and 8 diameters downstream was 
obtained for Mach numbers of 0.3 and 0.6. We first present an overview of the theory 
involved. We finish up with a statement of the work performed and finally provide 
charts from a 1999 APS talk which summarizes the results. 


2 Introduction 

The complementary technique consists of projecting the POD eigenfunctions onto 
an estimated velocity field obtained from application of LSE as described by Cole et 
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al [2] to obtain estimated random coefficients. These estimated random coefficients are 
then used in conjunction with the first couple of POD eigenfunctions to reconstruct 
the estimated random field, i.e., a low dimensional description. Bonnet et al 1994 
presented a qualitative comparison between, the first POD mode representation of the 
estimated random velocity field, and that obtained utilizing the original measured 
field, and found that the two are remarkably similar, in both the axisymmetric jet 
and two dimensional mixing layer. In order to quantitatively assess the technique, the 
root mean square (RMS) velocities were computed from the estimated and original 
velocity fields and comparisons made. In both flows the RMS velocities captured 
using the first POD mode of the estimated field are very close to those obtained 
from the first POD mode of the unestimated original field. These results show r ed that 
the complementary technique, which combines LSE and POD, allows one to obtain 
time dependent information from the POD while greatly reducing the amount of 
instantaneous data required. Hence, it is not necessary to measure the instantaneous 
velocity field at all points in space simultaneously to obtain the phase of the structures, 
but only at a few select spatial positions. Hence this type of an approach can be used 
to obtain experimental low dimensional descriptions of the jet. The w r ork of Bonnet et 
al 1994 applied the estimation in the radial direction only. In this study the method 
was extended to include the azimuthal direction. A brief review of the POD and LSE 
are included below followed by an overview of the complementary technique. 

2.1 POD Theory 

The POD was proposed by Lumley in 1967 as a mathematically unbiased technique 
for examining coherent structures in turbulent flows. He proposed that the coherent 
structure is the structure which has the largest mean square projection on the velocity 
field. In the following equations x denotes r and dx denotes rdr for the jet. If <f>(x, t ) 
is taken to be the candidate structure and then projected onto the velocity vector 
field, u(x,t), in the following manner, 

\u ■ (p\ 2 = \a\ 2 (1) 

a resulting structure which maximizes energy can be chosen. Equation 1 is assumed 
to be normalized by the modulus of $(x, t) since it is the degree of projection not the 
amplitude that is of interest in this study. Through the use of variational calculus 
this projection process can also be written as, 

j Rij(x, x 1 , t, i')^ n) (f' I t')d£dt' = \ {n) 4 n) (x, t ) (2) 

where the kernel of this integral eigenvalue problem is the cross-correlation tensor, 
Rij{x, x', t, t'). Since R X j is a symmetric function the solutions can be discussed in 
terms of the Hilbert-Schmidt theory(v. Lumley[4]). This implies that there is not one, 
but an infinite number of orthogonal solutions(eigenmodes) which can be normalized 
such that, 

J 4 n) ( x, t)</>S m) (:r, t)dx = 5 nm . (3) 
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( 4 ) 


The original random vector field can be reconstructed in terms of <f>\ n '* as follows: 

OC 

Ui{x,t) = 

71=1 

where the random coefficients a n (t ) can be calculated from: 

a n (t) = j u l (x,t)(f>\ n) {x,t)dx. (5) 

The contribution from each eigenmode to the turbulent stresses can be determined 
as follows: 

uiuj = \ {n) (f>i n) (x, t). (6) 

71 

In the jet study the spectral tensor may be defined by the following equation, 

S„ (r,r‘, f, m,z) = j Kj (r, /, r, 9, (7) 

where / denotes frequency, m denotes azimuthal mode number, r is the separation in 
time, 9 the separation in the azimuthal direction and z represent streamwise locations 
where the correlation tensors were measured (in the Ukeiley data, 4, 5, 6, 7, 8 diameters 
downstream for Mach numbers of 0.3, 0.6 and 0.85). Equation 2 now becomes 

J Sij(r,r\ f,m,z)r'dr' - A (n) (/,m)^- n) (r>/) m .^)» ( 8 ) 

where the ip's are frequency and azimuthal mode number dependent eigenfunctions 
and now represents the eigenspectra. 

The numerical approximation, detailed by Glauser et al[3], simply consists of 
replacing the integral in equation 8 by an appropriate quadrature rule (in this study a 
trapezoidal rule). 5 t j(r, r', /, m, z) is then obtained from experimental measurements 
of Ukeiley for the various downstream locations and utilized in equation 8 to obtain 
the eigenvalues and eigenfunctions. These eigenfunctions and eigenvalues will then 
be used in the complementary technique as described below. 


2.2 Stochastic Estimation Theory 

In the following equations x denotes r in the jet. Only estimates in r are shown 
here. It is straight forward to extend the analysis to include additional directions as 
shown by Cole et al (see Cole et al APS 1997 and www.clarkson.edu/glauser) in their 
application to the time dependent PIV (for this study it will be extended to include 
the 9 direction). In general, the conditional average is defined as 

g(u)\E — expected value of g(u ) (9) 

provided that the event E , the detector of coherent structure, occurs. The properties 
of these coherent structures are not known beforehand, therefore it is difficult to select 
reliable unambiguous and unbiased detector events. Adrian [1] suggested choosing 

g(u) = u(x') (10) 
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and 


E = e < u(x ) < e + de. (11) 

This confines the velocity vector to a small interval between e and e + de where e is 
any arbitrary vector. This can be expressed as 

u{x') = u{x')\u{x) (12) 

which can be approximated, using a Taylor series expansion about u{x) = 0, as 


Ui(x') = Aij(x')uj(x) + Bij k (x')uj(x)u k (x) + Cijki{x')uj(x)uk(x)ui{x) H (13) 


Values for the estimation coefficients, A lJ (x'), Bij k (x'), Cij k i(x') are selected such that 
the mean square error is minimized, ie., 


e { = [ui(x') - Ui(x')\u(x )] 2 


for i = 1, 2, 3. This minimization requires that 

dej _ dej _ dej _ 

dAijlx') dB ljk (x') dCijki(x') 

which produces the following equation 

u J (x)u k (x)A i k{x') + Uj(x)u k {x)ui{x)B M {x') 

+U j (x)u k {x)u l (x)u m (x)Ciklm(x') = Uj(x)Ui(x'). 


(14) 


(15) 


(16) 


2.2.1 Linear Stochastic Estimation 

Tung and Adrian [5] have shown that linear stochastic estimation produces reasonable 
qualitative estimates and little is to be gained by using second order or higher. By 
applying linear stochastic estimation only the first term on the right hand side of 
equation 13 is retained. As a result, equation 13 becomes 

Ui(x') = Aij(x')uj(x ) (17) 

This also reduces equation 16 to: 


Uj(x)u k (x)A ik (x') = Uj(x)ui(x') (18) 

where Uj(x)u k (x) is the Reynolds stress tensor and Uj(x)ui(x') is the two-point cor- 
relation tensor. 

Linear stochastic estimation is a useful tool for the identification of coherent struc- 
tures in the axisymmetric jet mixing layer as has been demonstrated by Cole et al [2]. 
After applying this technique to u,v data, the matrices that are a direct result of the 
expansion of equation 13 for a two wire estimate are: 
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First System: 



Vrefi Vrefi 

V 2 

V refi 

Vrefi V r e / 2 

U re f l V re f 2 


r 4 re /i 1 

^ 1 1 u; 

A refl 

^12 w 
4 re fo 


Vrefi Vw 

V re fi ^re/i 

Vrefi Vre / 2 
U ref2 

Vrefi Vref2 


— 

Vrefi Vw 

lire fi^re f\ 

V're f^Vre }\ 

Vref2 Vre / 2 


A 

J re /2 

L ^12 w J 


Vref2Vw 




y 2 . 
ref 2 J 



_ V re f 2 V w _ 

_ V re f 2 U r efi 

Vref^Vre f\ 

V re / 2 V re f 2 




Second System: 

U re.fi ^refi Vref\ ^ re /j tfre/2 

^re/i^re/i ^re/i ^re/i^re/2 

U r ef 2 ^refi V re j 2 V re j x U ref 2 

^re/ 2 ^re/i V re j 2 V re f x V re f 2 U re f 2 


r 

5- 

53 


r A refi 1 

^2\w 
A re fl 
^22 w 

\ re h 

^2lw 

A ref 2 

L ^22u; J 


Vrefi Vw 

Vrefi Vref 2 


— 

Vrefi Vw 

Vre / 2 Vre / 2 



Vre / 2 Vyj 

V ref 2 



_ Vref2 Vw _ 


( 19 ) 


( 20 ) 


where re/i and re/ 2 refer to reference wires 1 and 2 respectively, and u; refers to 
the wire number. It should be noted, that for these systems of equations, only the 
two-point space-time correlation data is utilized. These systems are not a function 
of the condition being investigated. The estimated velocity components for the two 
wire reference case can then be found from the expansion of equation 17 ie. , 

U v) = A\\^uc re{l + A\%vc refl + A\\ f *uc Tef2 + A r y^ vc reh (21) 

and 

v w — ''f 21 u! uc ref\ + A T ^/w VC Te f l + A r 2 \^UC rc f 2 + A 2 2wVC re f 2 . (22) 


It is in these estimated velocity equations that the condition selected plays a role 
(i.e., through uc re j 1 , uc re f 2 , vc Te f l and vc re j 2 ). A single wire estimate is obtained by 
merely setting all terms containing re/ 2 = 0. Without much trouble this system can 
easily be expanded to include addition wires distributed in the 9 direction as done in 
this study. 


2.3 Complementary Technique 

Mathematically the stochastic estimates of the random coefficients are calculated for 
a given location downstream from: 

a e n st (f,m) = J u* s \rj,m)i>\ n) \rj,m)rdr (23) 

where u^ st (x, f,m) is either a single or multipoint linear stochastic estimate of the 
random field and i />{ n ^*(r, /, m) is obtained from the original POD eigenvalue problem. 
The estimated streamwise or radial velocity can be reproduced in Fourier space by 

OO 

u e i st (x, f,m) = a e n st (f, m)V>i n) (r, /, m) (24) 

71=1 

and then inverse transformed to obtain u® st (r, 0, f). 

In this study we have been concentrating on the first POD mode representation 
u\^ est (r,9,t), the low dimensional discription. 
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3 Work Performed 


3.1 Experimental Low Dimensional Description via LSE 

We ran experiments at the downstream locations where the two point hot wire mea- 
surements were obtained by Ukeiley and applied the LSE. These experiments involved 
sampling simultaneously hot wire signals from a relatively course spatial grid in r 
and 9. From this simultaneous data, coupled with the two-point measurements of 
Ukeiley via the LSE, an experimental low dimensional description of the jet at 4, 5, 
6, 7 and 8 diameters downstream was obtained for Mach numbers of 0.3 and 0.6. 
An APS talk is attached which summarizes the work. Also reference the web site 
www.clarkson.edu/ etrl/taylorja/CompressibleJet/index.html for more details and 
animations of the LSE time series. 
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Introduction - Motivation 


SUB-SONIC JET 
NOISE REDUCTION 


*** Requires *** 

An Improved Understanding of the 
Turbulent Jet Mixing Layer 


Introduction - Some Previous Studies 


Arndt, R., George, W.K. 

What's the relationship between Coherent 
Structures and Jet Noise? (1972) 


Amdt, R., and Long, D. - University of Minnesota 


George, W.K. and Glauser, M.N.- SUNY at Buffalo/Clarkson University 
x/D: 3 

Re # - 100,000 
Ma # ~ 0.0 

13 single component hot-wire or 4 x-wire probes 

George, W.K. and Citriniti, J.H. - SUNY at Buffalo 

x/D: 3 

Re # ~ 80,000 
Ma # ~ 0.0 

138 single component hot-wire probes 

Ukeiley, L.S. and Seiner, J.M. - NASA Langley/University of Mississippi 

x/D: 4,5, 6, 7, 8, 9, 10, 11, 12 
Re # ~ 300,000; 600,000; 1,000,000 
Ma # ~ 0.3, 0.6, 0.85 
12 X-Wire probes 

Current Study - Clarkson University/University of Mississippi 
x/D: 4, 4.5, 5, 6, 7, 8 
Re # ~ 300,000; 600,000 
Ma # ~ 0.3, 0.6 
15 X-Wire probes 


Linear Stochastic Estimation - Brief Review 


Adrian, R.J. (1977) 

Proposed the use of “Stochastic Estimation” as a 
means of using conditional information on a course 
grid of points to estimate, or infer, the behavior of the 
flow on a finer grid. 


Cole, D.R., Glauser,M.N., and Guezennec,Y.G. (1992) 

Proposed using the instantaneous fluctuating velocity 
on a course grid as the input variable, or condition, in 
the estimation as a means of reducing experimental 
complexity. 

Bonnet, J.P., Cole, D.R., Deville, J., Glauser, M.N., 
and Ukeiley, L.S. (1994) 

Proposed the “Complementary Technique” as a 
combination of the LSE and POD. 

1 . Obtain two-point correlation tensor 

2. Solve the POD for the eigenfunctions/eigenvalues 

3. Use the LSE to obtain an estimate of the 
instantaneous flow field 

4. Project the eigenfunctions onto the flow field to 
generate an experimental low-dimensional 
description of the flow 


Experimental Description - Facility 


Small Anechoic Jet Facility (SAJF) 


Features: 

Continuous Operation at 2 lbm/sec (0.9 kg/sec) 
Electric Heat up to 1000 °F (538 °C) 

Exit Turbulence Intensities « 0.15% 



Coflow ~ 2% of main flow 




Experimental Description - Two-Rake Experiment 


A 


r/D = 2.00 

1 


1 


i 



Experimental Parameters: 

2 Rakes of Probes 

12 Auspex X-Wire probes 
19 Azimuthal Locations (0°:10°:180°) 

6 Radial Locations (r/D: 0.175, 0.339, 0.504, 0.668, 0.833, 0.999) 

3 Mach Numbers (0.3, 0.6, and 0.85) 

Data Acquisition Parameters: 

65 kHz Sampling Frequency 

31.5 kHz Low-Pass Filter Cut-Off Frequency 

16 bit A/D Conversion 

256 blocks with 2048 samples/block 





Experimental Description - LSE Experiment 



Experimental Parameters: 

Single Probe Holder on a Course Grid 

15 Auspex X- Wires Probes 

12 Azimuthal Locations (0°: 30°: 360°) 

2 Radial Locations (r/D: 0.175, 0.504) 

Data Acquisition Parameters: 

65 kHz Sampling Frequency 

25 kHz Low-Pass Filter Cut-Off Frequency 

16 bit A/D Conversion 

100 blocks with 2048 samples/block 


v-component of mass-flux kg/(m s) u-component of mass-flux kg/(m s) 


Results - Time Histories: Mach = 0.3, x/D = 4 




Results - Time Histories: Mach = 0.6, x/D = 4 
















Results - Movies 


Citriniti, J.H.: Instantaneous Velocity 
Ma # -0.0 
Re # -80,000 
x/D: 3 


Current Study: Estimated Mass-Flux 
Ma # -0.3 
Re # -300,000 
x/D: 4 



Current Study: Estimated Mass-Flux 
Ma # -0.6 
Re # -600,000 
x/D: 4 



Future Work - 

Experimental Low-Dimensional Description 
using the Complementary Technique 


Citriniti, J.H.: 

1 POD Mode; 0,3, 4, 5, 6 Az. Modes 
Ma # -0.6 
Re # -600,000 
x/D: 4 




Summary 


Jet mixing layer has a very similar multi-point, statistical 
behavior at: 

Ma # ~ 0.0 Re # -100,000 
Ma # ~ 0.3 Re #-300,000 
Ma # - 0.6 Re #-600,000 

(i.e., Azimuthal integral length scales are similar ) 


Time dependent behavior also exhibits similar behavior. 


(i.e., coherent in the potential region, with a higher 
azimuthal mode structure towards the outside of the 

shear layer) 


Future Work 


How does the l ow-dimensional description of the axi- symmetric 
shear layer change with Reynold’s number? Mach number? 


